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byline: Soliton-type Solutions and A^-component Systems 

Abstract 

The algebraic geometric approach to A^-component systems of nonlinear integrable PDE's is 
used to obtain and analyze explicit solutions of the coupled KdV and Dym equations. Detailed 
analysis of soliton fission, kink to anti-kink transitions and multi-peaked soliton solutions is 
carried out. Transformations are used to connect these solutions to several other equations that 
model physical phenomena in fluid dynamics and nonlinear optics. 

1 Introduction 

The solution of nonlinear evolution equations using techniques from algebraic geometry was ini- 
tially developed to handle A^-phase wave trains. With this approach, parameterized families of 
quasi-periodic and soliton solutions are associated with Hamiltonian flows on level sets of finite- 
dimensional phase spaces. In Section ||, these flows are described using the /i variable representation 
on symmetric products of Riemann surfaces. The first integrals in the quasi-periodic case have the 
form = C{fij) where C{^j) is a polynomial with constant coefficients and Pj is the conjugate 
variable for fij. The polynomial C{E) is called the spectral polynomial and determines the form of 
the first integrals. 

The algebraic geometric approach provides a way to construct solutions by analytical or nu- 
merical integration of a system of Hamiltonian equations for these /i variables. Solutions of the 

'PACS numbers 03.40. Gc, ll.lO.Ef, 68.10.-m, AMS Subject Classification 58F07, 70H99, 76B15 
■I^Research partially supported by NSF grant DMS 9626672. 
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nonlinear PDE's are expressed in terms of these variables by using the trace formulas. The fj, vari- 
able representation yields an action angle representation on a Jacobi variety (invariant variety in 
the phase space) that linearizes the Hamiltonian flow, and in the quasi-periodic case the solution 
of the /X equations is reduced to a Jacobi inversion problem. Solutions are then expressed using 
Riemann theta functions. For details see, for example, Mumford(1983) and Ercolani and McKean 
(1990). 

As pairs of roots of the polynomial C{E) coalesce, soliton solutions begin to appear. Applying 
this limit to the first integrals and the equations of motion in terms of the jj, variables for quasi- 
periodic solutions, Hamiltonian systems of ODE's that describe soliton solutions are obtained. 
Soliton solutions are computed by solving these equations and using the trace formula to connect 
them to the associated nonlinear PDE's. In the soliton limit the Jacobi inversion problem for 
the system often reduces to a system of algebraic equations. As will be shown, these algebraic 
equations are exactly solvable in the case of genus one and two. For details about the connection 
between soliton and quasi-periodic solutions see for example Ablowitz and Ma (1981) and Alber and 
Alber(1985). The exact representations obtained for soliton solutions are shown to be related to 
the Hirota r-functions. The case of genus n solutions is analogous to the lower genus case and can 
be described using similar formulas. Using the action-angle representation, the algebraic geometric 
approach also introduces a powerful way to compute the phase shifts due to soliton interactions 
(see Alber and Marsden(1992) for details.) 

In this paper, soliton solutions of multicomponent systems of equations are studied using the 
algebraic geometric approach. The soliton fission effect, kink to anti-kink transitions, and multi- 
peaked solitons are demonstrated using a class of commuting Hamiltonian systems on Riemann 
surfaces. These first two effects manifest themselves in the soliton limit of the genus two quasi- 
periodic solution when six roots of the spectral polynomial C{E) coalesce in a pairwise fashion. 
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Exact formulas for these solutions are obtained and asymptotic and numerical analysis of them 
is performed. The technique used to obtain these limiting solutions is demonstrated explicitly in 
the case of the coupled Korteweg-de Vries (cKdV) and the coupled Dym (cDym) equations and 
physically relevant equations associated with them. 
The modified cKdV system is 

Ut = Vx - ^UUx + KiUx , (1.1) 
Vt = \uxxx - VUx - \uvx + KiVx , (1.2) 

which reduces to the cKdV system for Ki = 0. Otherwise Ki is viewed as a small constant. The 
coupled Dym equations are given by 

Ut = \uxxx - ^UUx + Vx + KiUx , (1.3) 
Vt = -UxV - \uvx + KiVx . (1.4) 

The general method is demonstrated by describing quasi-periodic and soliton solutions for the 
cKdV and cDym systems for genus two and less. 

The cKdV and cDym equations are both generic examples of A^-component systems. Energy 
dependent Schrodinger operators and 6z-Hamiltonian structures for multicomponent systems were 
investigated in Antonowicz and Fordy(1987). Quasi-periodic and soliton solutions were studied in 
connection with Hamiltonian systems on Riemann surfaces in Alber et a/. (1997). In Alber et al. 
(1994) it was also shown that the presence of a pole in the associated Schrodinger operator yields 
a special class of weak billiard solutions for nonlinear PDE's. 

The soliton fission effect, kink to anti-kink transitions, and multi-peaked solitons extend to 
equations that model physical phenomena. The generalized Kaup equation, the classical Boussinesq 
system, and the equations governing second harmonic generation (SHG) are each connected to the 
cKdV system through nonsingular transformations. Direct application of these transformations 
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enables solutions of the cKdV system to be interpreted in the context of these related equations. 
Such transformations are given explicitly in Appendix Both the Boussinesq System, 

Ut + + UU^ = -fU, , (1.5) 
Wt + U,., + (WU), = ^W, , (1.6) 

and the generalized Kaup equations 

n = <l>xx + li'i^ - ^(^)S'^(pxxxx - e{(t)xTr)x + aiTx , (1.7) 
(f)t = -acl)x + ^4)1 + TT , (1.8) 

arise from the theory of shallow water waves (see Whitham(1974)). Here 7 and a are small param- 
eters. 

In optics, the interaction of a wave envelope at frequency and wavenumber {w, k) with a second 
wave at twice the frequency is modeled by the system of equations 

{qi)x = -2(?29i* , (1-9) 
{q2)r = ql . (1.10) 

This process is called second harmonic generation in nonlinear optics and is used to convert laser 
light to its second harmonic frequency. 

The scattering problem for the energy dependent Schodinger operators was studied by Jaulent(1972) 



and Jaulent and Jean(1976). The completely integrable variant of the Boussinesq system (1.5)- 
(|l.6| ) was first introduced by Kaup(1972). In Matveev and Yavor(1979) ^-functions were used 
to describe quasi-periodic solutions of the Boussinesq System. They also described a particular 
type of A^-soliton solution using singular classes of 0-functions. Rational solutions were studied 
in Sachs(1998) in connection with a pair of Calogero-Moser equations coupled through the con- 
straints. Martinez Alonso and Medina Reus(1992) and Estevez et a/. (1994) described some of the 



soliton solutions using Hirota's r-functions. Using asymptotics of these soliton solutions they also 
demonstrated soliton fission. A connection between the SHG system and the cKdV system was 
recently discussed by Khusnutdinova and Steudel(1998). 



2 Generating Equations for the Coupled KdV and Dym Equations 

We begin by describing the general approach of generating equations and applying it to the cKdV 
and cDym equations. Details for the general case of A^-component systems are discussed in Alber 
et a/.(1997). 

2.1 Dynamical Generating Equations. 

The hierarchy of the coupled KdV and Dym equations is obtained as the compatibility condition 
for the eigenfunction of the linear system of equations 



The time flow is produced by the linear differential operator 

where B{x,t,E) is a specified rational function. The operator L is assumed to be of the energy 
dependent Schrodinger type, 



LtP = , 



(2.1) 



ipt = Alp . 



(2.2) 



(2.4) 



with a rational potential having the form 



V{x,t,E) 




(2.5) 
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where rj are constants and Vj{x,t) are functions of the variable x and the parameter t. E is a 
complex spectral parameter. In particular, the potential is chosen as 

V{E) = kE^ + u{x,t)E + v{x,t) , (2.6) 

for the cKdV system or 

V{E)=uix,t)+E + ^^ , (2.7) 

to recover the cDym system. Here k = ±1. One chooses n = —1 to establish the transformation 
from the cKdV system to the SHG system, and k = 1 to establish the transformation from the 
cKdV system to the Boussinesq System. Notice that the main difference between the cKdV and 
cDym systems is the presence of a pole in the Schrodinger operator ( |2.4| ) associated with the cDym 
equations. The pole in the potential for the cDym case was shown in Alber et a/. (1994) to be a 
necessary feature for systems with weak billiard solutions. 

The compatibility condition of (|2.1| )- (|2.2|) can be found by taking the t derivative of (0), acting 
on (|2.2| ) by L, and forcing the fact that these two operators commute. This leads to the following 
system of equations 

iLt + [L,A])i; = 0, Lij = 0, (2.8) 

where [L, A] = LA — AL is the commutator of L and A. Using the definition of the differential 

operator A in ( |2.3D and L in ( p.4D , this Lax equation yields 

dV _ Id^B ^dB^ ^dV 
dt 2 dx^ dx dx ' 

which is a generating equation for the coefficients of the differential operator A. By taking B to be 

the rational function, 

m n 

Bix, t,E)=Y^ bm-kix, t) E^ = E-^ \{{E- fikix, t)) , (2.10) 



substituting it into the generating equation (^^) and equating like powers of E, a recurrence chain 
of equations for the coefficients bj is obtained. Evaluating these equations one by one, a PDE for 
the coefficient 5„ is obtained, where n = m + r. By considering all possible values of n and m, a 
hierarchy of systems generated by the Lax equation with a given potential (generating equation) is 
obtained. 

Assuming that — — = in ( |2.9| ) and integrating gives the stationary generating equation which 
at — 



has the form 



-B"B + UB'f + 2B^V = C{E) , (2.11) 



where the choice of B{E) from ( 2.10| ) ensures that C{E) is a rational function with constant 



coefficients. These coefficients are the ffist integrals and parameters of the coupled system of 
equations. (For details about the general method see Alber and Alber(1985).) This equation gives 
rise to level sets in the phase space C" corresponding to the Riemann surface 

= C{E) . (2.12) 

The dynamics of genus n quasi-periodic solutions for u and v with respect to the x and t coordinates 
are captured as flows on the level set produced by a symmetric product of n copies of the Riemann 
surface ( |2.12| ) (see Alber et aZ.(1997) for details). 

The method of generating equations yields the cKdV and cDym equations by using B{x,t,E) = 
60(2;, t)E + bi{x, t). In Appendix |X| solutions of the generalized Boussinesq and generalized Kaup 
equations are linked to these systems. Further, the second harmonic generation equations are 
obtained when B[x,t,E) = b2{x,t)E^^ . 
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3 Solutions of the cKdV and cDym Systems from Dynamical Sys- 
tems on Riemann Surfaces 

In this section we obtain finite-dimensional Hamiltonian systems on Riemann surfaces for the /i 
variables defined in ( 2.1C| ) as the roots of the function B{x, t, E). These systems capture the essential 



dynamics of quasi-periodic and soliton solutions of the cKdV and cDym systems. The quasi-periodic 
solutions are often called n-gap solutions in physics literature and quasi-periodic solutions of genus 
n in mathematics literature. In Appendix ^ the solutions of these finite-dimensional Hamiltonian 
systems are linked to the solutions of the integrable nonlinear PDE's through trace formulas of 
the general form u = a Yl^=i + where a and /? are constants. This general construction also 
provides links between the solutions of the cKdV and cDym systems and other nonlinear PDE's. 
Choose B{x,t,E) as in ( 2.10| ) where n = r + m is the genus of the desired solution. For the 



cKdV case, we substitute ( p.6| ) into ( p.9[ ) and ( 2.11| ) to obtain 



utE + vt = -\B"' + 2kB'E'^ + 2B'uE + 2B'v + Bu' E + Bv' , (3.1) 



and 



-B"B + \{B'f + 2kB'^E'^ + 2B'^uE + 25^ = C{E) . (3.2) 
For the cDym system we substitute the potential ( |2.7| ) into the same equations to obtain 

utE + vt = -\B"'E + 2B'Eu + 2B'E'^ + 2B'v + Bu E + Bv' , (3.3) 



and 



2B'^v 

~e' 



-B"B + l{B'f + 2B^u + 2B'^E + = C{E) . (3.4) 



By equating like powers of E on the left and right hand side of the equations ( |3.2D and (3.4) we 



obtain the necessary forms for C{E). Therefore we write 



2{n+l) 

C{E) = 2kE-^'- n ~ "^^) ' (^•^) 
1=1 



for cKdV equations and 



2{n+l) 

C{E) = 2^-(2-+i) n - "^^) ' (3.6) 
1=1 

for the cDym equations, where the real numbers mj are the roots of the polynomials C{E). 
3.1 Finite-Dimensional Hamiltonian Systems. 

From the generating equation above it follows that the fii's, which are the roots of the function 
B{x,t,E), are solutions of finite-dimensional Hamiltonian systems. By solving these Hamiltonian 
systems and using the trace formula, the dynamics of the roots fii are connected to the functions 
u and V to obtain solutions of the PDE's. 

Namely, Hamiltonian equations for /Xj's are obtained by substituting E = fii into the generating 
equations ( |3.2D , ( 13.4]) . This yields the following systems of equations for cKdV: 



w ^lj=l ""J J 

lii=± ^ ^ z = l,...,n, 3.7) 



for the flow in space, and 



ii — rrij 



= =F T=f 7 ^ i = l,...,n, (3.8) 

for the flow in time. Here n = m + r from the definition of B in (|2.1C| ) and rrij are fixed real roots 



of C{E) in (3^). The plus/minus refers to which branch of the Riemann surface ( 2.12 ) the solution 



is on. Systems ( ^.7D -(p^) are commuting Hamiltonian systems with Hamiltonians 
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where D[^i) = 1 and D[fii) = — Yljjki f^j the stationary and dynamical cases respectively. These 
systems share the same complete set of first integrals: Pj = C{fj,j),j = 1, n where the polynomial 
C{E) is defined by Q (for details see Alber et a/.(1997)). 
For the cDym equation we see that 



fii = ±^—^ r z = l,...,n, (3.10) 



for the spatial flow, and 



= =F ^=f^ 7 ^ i = l,...,n, (3.11) 



for the time flow. Systems ( 3.10D -( ^.ll ) are also Hamiltonian systems and they share the same 



complete set of first integrals: = C{fij),j = 1, ...,n where polynomial C{E) is defined by ( |3.6| ). 
Systems (^)-([3^) and ( 3.10[ )-( pT ) can be solved analytically by reducing them to Jacobi 



inversion problems. We will demonstrate the general method in the next section. These equations 
are also easily integrated numerically. Perhaps the best method for accomplishing numerical inte- 
gration with the use of symplectic integrators. Such integrators preserve the Poincare invariants 
and are stable over a long period of time, see for example Channell and Scovel(1990). 

3.2 The Trace Formulas. 

The connection between the solutions u and v of the cKdV equation and the /Hj's from (l3.7|) ,(^) 
is derived in Appendix ^ and is given by 



u 

1=1 



2K^lii + 2KKi , (3.12) 



-2k fiifij + - Km + K2 , (3.13) 



4fi; 

i<j<n 
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where 

^ 2n+2 

Ki = "2 E ' (3.14) 

i=l 

K2 = K mirrij — kK^ . (3.15) 

i<j<2n+2 

Observe that the smah parameter Ki from ( |1.1| ) is zero if 'Y^=i'^ ~ ^- '^^^ trace formulas for 
the cDym system are 



u 

i=l 



2^^i + 2Ki, (3.16) 

i=i 

V = -^u" - 2 fiifij + ^u^ - Kiu + K2 . (3.17) 

l<i<j<ra 

4 Classification of Limits of Periodic Solutions 

For the next three sections we seek various periodic (genus-one) solutions of the Boussinesq System. 
Therefore unless otherwise stated we assume r = from (|1|) and At = 1 from (U). We first 
obtain the periodic traveling-wave solutions and then show that they are equivalent to solutions 
obtained in terms of a variable in the case n = 1. This provides a natural introduction to the 
algebraic geometric method. A one-soliton solution is then obtained by deforming the Riemann 
surface of the genus-one periodic solution. This method of first finding periodic/quasi-periodic 
solutions and then deforming the level set (Riemann surface) in the phase space to obtain soliton 
solutions will be utilized throughout this paper. (For details about general approach see amongst 
others Ablowitz and Ma(1981) and Alber and Alber(1985).) 
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4.1 Periodic Traveling- Wave Solutions 

Let U = U{C) and W = W{C) where C = 2; — ct so that the Boussinesq System becomes 

-cU' + W' + UU' = jU' , (4.1) 

-cW + U'" + {WUy = -iW' , (4.2) 

where W and U' denote differentiation with respect to Then (|4.1|) gives 



W = r]U -W'^ + T0 , (4.3) 



where 77 = c + 7 and tq is a constant of integration. Plugging this into ( f4.2D and integrating twice 
gives 

^-^U' + + i([/')2 - = nU + r2 , (4.4) 

where ti,T2 are constants of integration. Writing this as an integral equation and taking a square 
root we obtain 

dC = ± (4.5) 

where 

4 

C4(C/) = Y[{U - rrij) = U^- ArjU^ - 4(ro - t]'^)U^ - SnU - Sra- 
1=1 

Notice that the right hand side of this differential equation is multi valued since it involves a square 
root. This is uniquely defined on a Riemann surface of genus 1 parametrized by a pair (W, E) 
where 

4 

W'^ = Ci{E) =J\{E -mj). (4.6) 
1=1 

One indicates one of two sheets of the Riemann surface by choosing a particular sign in front of 
the square root: W = ±^Ca{E). Therefore U in ( |4.5| ) is considered on a particular sheet of the 
Riemann surface ( [4.6| ) and so meaningful integration can take place. 
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The equation ( [4.5[ ) can also be obtained from the /U-equations ( |3.7D and ( |3.8[ ) for n = 1. Notice 
that the trace formula shows that U and /i are linearly related so that one can substitute ^ instead 
of U in ([4.5|). Using /i in this case will make it consistent with the formulas in the case when n = 2 
to be described in the next section. 



After integrating (4^), the following Jacobi inversion problem is obtained 



e = x-ct + eo = \ r . = , (4.7) 

As stated before, this is a typical Jacobi inversion problem (Mumford (1983)). This integral is 
inverted using Jacobi's elliptic functions, 



(mi — m^){m2 — m.4) dfj, dz 

2 7/^0 



(4.5 



where 



2 _ ("^2 - "i4)(^ - mi) 2 _ ("^1 - rni){m2 - 1713) 



(mi — m4)(;U — m2) ' (m2 — m4)(mi — ms) 

Notice that ( |4.8| ) is an elliptic integral of the first kind (see Mumford (1983)). Therefore, n is 
obtained explicitly and 

_ m2(m4 - mi)sn^(fc, uj) + mi(m2 - m4) 
m2 — m4 + (m4 — mijsn^(/e, to) 

where sn(A;,a;) is the Jacobi sine function and oj = {x — ct + Oq) j ^ (mi — m3)(m2 — m4). This 



function is plotted in Figure B.l 



Periodic solutions correspond to the case when all m^ are distinct. In this case, each fixed 
point in the phase space repels the trajectories so that for any initial value, /i oscillates periodically 
between the two nearest points. Notice that /i(0) must be chosen so that the right hand side of ( [4. 71) 
is real. Then the fixed points repel and [i remains real valued. The solution leads to solutions 
for \J and W through the trace formulas. The shape of U is essentially the same as /i. The shape 
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of W is aperiodic and is discussed further in the next section. Equation ([4.7|) can be interpreted as 



defining an angle variable 9 where c is then the action variable. From (4.7) it follows that in terms 
of these variables the initial Hamiltonian flow linearizes. It also can be viewed as an Abel -Jacobi 
map from a hyper elliptic curve 

= C{E) = 2{fi - mi)ifi - m2){li - m^){ii - m^) , (4.11) 

or in other words, a Riemann surface of genus one, onto the Jacobi variety: J = [C|tt;Z] where C 
is a complex plane and w is the period lattice of the holomorphic differential from ( |4.7[) . 

4.2 One-Soliton Solution of Kaup Type 

To examine soliton solutions, one deforms the Riemann surface W"^ = Ci{E). (For details about 
soliton deformations see amongst others Ablowitz and Ma(1981) and Alber and Alber(1985).) Here 
we consider the limit mi — > m2 — > a, i.e. where the 2 roots mi and m2 coalesce into one point. 
As this limit is approached, a soliton solution is obtained as the period of the periodic solution 
increases to infinity. In this case a is called a double point. This double point is an attractor in the 
phase space. Without loss of generality, we assume a < m^ < 7714. This assumption leads to the 
realization that the only pertinent solution is obtained when a < /i(0) < m^. 
On the principle branch of the square root, 

x-ct + Oo = I r = , (4.12) 

= - [(a - m4)(a - m^)]-'/' arctanh f j i° ~ "^'^^ ~ , (4.13) 



so that 



(a - m'i){n - mi) ) 



^ _ msjm^ - g) + m^ia - m3)tanh^(a;) 

(777-4 — a) + {a — 7n,3)tanh^(u;) ' 
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where uj = —{x — ct + 9o)\/ {a — 77^4) (a — ms) and a, m3,m4 are functions of c. Using the trace 
formulas ( p.l2D -( |3.13D and the transformations ( |A.3| )-( |A^ ), the exact formulas for U and W from 
the Boussinesq equation are obtained. For example, 

_ -^a? + (ms - 771,4)^ + 2a(m3 + 771,4) + - m3)cosh(2u;) 
2a — (m3 + 7714) + (?77-3 — 777,4) cosh (2cj) 

Notice that /x has a shape similar to a KdV solitary wave. It has a peak at 7773 and approaches 
7771 for large |rE|. From this we conclude that U is also shaped like a KdV soliton with a peak of 
height 4a + 27774 — 27773 and for large \x\ approaches 27713 + 27774. One might be surprised that the 
soliton does not approach as x — > =boo. But remember that this is the modified cKdV equation 
with Ki 7^ 0. Howver we may choose the parameters so that the solution does approach zero. The 



solution is plotted in Figure B.2. This solution was first found by Kaup(1972) using the inverse 
scattering transform in case when 7773 = —777,4 and a = 0, hence Ki = 0. One advantage of the 
approach used here is that W is easily found using the trace formulas, and the relationship between 
U and W is seen explicitly. 

4.3 Solutions with Two Peaks. 

Observe next that since n = 1, is only a quadratic in U. The polynomial 

W = -^U"^ + (2a + 7773 + mi)U - iarris - ianii + (7773 - 7774)^ , (4.16) 

obtained from the trace formula is a parabola in U with vertex at U = 2a + 7773 + 7774. Since U has 
the shape of a solitary wave, W has either one or two peaks depending on the parameters defining 
the vertex of U. (See Figure |B.2| ). If U is concave down, then is a double peaked soliton if 
7773 + 7774 < 2a and 87773 < 2a + 7774. In the unperturbed case (7 = 0), these conditions reduce to 
7773 < < a. Otherwise W has a single peak. The concave up case is similar. 
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4.4 One Kink Solution 

The one-kink solutions were initially found by Alonso and Rues(1992) as the simplest solutions 
obtained using the bilinear formalism of the Kyoto school. Our method introduces an alternate 
description of soliton fusion and fission for n = 2 and simplifies the computations. 

To obtain the one- kink solutions, we deform the Riemann surface further by taking the limit 
ms m4, so that rrii = m2 = ai and = = 02. Note that if k = — 1 this limit is not 
permitted for real valued fi in ( p. 7] ). However, for k = 1 the system is readily solved for real fi 
producing a special case of equation ( |4.5| ) where ti = = r2 = 7. After inverting, the angle variable 
is 

e = x-ct + eo = ±1 jjj^^ = T 'arctanh (^) , (4.17) 

and we find using the trace formulas and the transformations obtained in Appendix |^ that 

U = ±ctanh(-|(x -ct + 6lo)) + C2 , (4.18) 
2 

W = ±jsech^{-'^{x-ct + 9o)) . (4.19) 



These functions are plotted in Figure B.2. Here U is always a kink or anti-kink while W is always 
similar to a typical KdV soliton. Observe that the speed of the soliton, is exactly the amplitude 
of W and is proportional to the square of the amplitude of the U soliton. This connection between 
wave speed and amplitude is reminiscent of that found in the KdV equation. 

5 Classification of Limits of Genus-Two Solutions 

In this section genus-two quasi-periodic solutions are constructed. By deforming the Riemann 
surface (spectral polynomial), two-soliton solutions are obtained and several types of soliton-soliton 
interactions are described. 
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5.1 Quasi-periodic Genus-Two Solutions 

In the genus-two case U is given by the trace formulas to be the sum of two periodic functions, 
Hi and ii2, and a constant. In general the functions, /xi and /i2, have noncommensurate periods, 
so that the solution U is generically quasi-periodic and systems (|3.7D and (^) are defined on the 
symmetric product of two copies of the Riemann surface (hyperelliptic curve) of genus two given 
by 

W^ = CeiE), (5.1) 

where 

6 

Ce{E) = kII{E - m,) . (5.2) 

1=1 

After reordering the equations, summing them up, and integrating, the following Jacobi inversion 
problem is obtained: 

" +r /^^ =2x-2a,t + gg, (5.3) 



^^^^^ + r =2x-2a,t + g°. (5.4) 



•4 \/CWM Jf^'^ \/C6(/i2) 

Inverting the Abel- Jacobi map defined by ( |5.3| )- (|5.4| ) results in expressions for fii and /i2 in terms 
of Riemann ^-functions. (For details about the Abel- Jacobi map see Mumford(1983), Matveev and 
Yavor(1979), and Ercolani and McKean (1990).) 

Having derived the genus- two quasi-periodic solutions, several limiting cases will now be ex- 
plored to introduce solitons. Below each distinct case is considered. 

5.2 One-Soliton Solution on a Quasi-periodic Background 

A one-soliton solution on a quasi-periodic background is obtained in the limit mi m2 a. Just 
as in the one-soliton case, a soliton is created as the Riemann surface is manipulated by pinching 
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two elements of the spectrum. However, in this case there are two fj, variables and the orbit for only 
one /i-variable is changed. The other fi variable remains periodic. The result is a solution with a 



Kaup type soliton on a quasi-periodic background, and it is plotted in Figure p. 4 The problem of 
inversion may be written in the following way, 

+ r = 2t + , (5.5) 

(/^i -a)\AfMm) J^^'^ (/i2 -a)x/^4(/i2) 

dill ff"^ dii2 , X 

^ + / ^ =2x- 2at + , (5.6) 

where 



Pi{E) = {E- m^){E - m^){E - m^){E - mg) . (5.7) 
(For details about inverting problems of this type see Alber and Fedorov(1999).) 

5.3 Two-Soliton Solutions of Kaup Type 

The two-soliton solution is obtained by piecewise pinching together two pairs of elements of the 
spectrum so that mi — > m2 ai and ms m/i a2- Here there are two double points ai 
and a2, as well as two remaining hyperelliptic points at and uiq. Choosing the initial data on 
the positive branch of the Riemann surface for both variables leads to the following problem of 



inversion 



Ml 



djJ-i 



{fii - a2)x/(/ii - m5)(/ii - me) 



d^l2 



+ 1 ^ = = 2x-2ait + g°, (5. 

4 (Ai2 - a2)\/(Ai2 - ?Ti5)(;U2 - me) 



'^i dm 



f,o (/ii - ai) Y^(/ii - m5)(//i - me) 

+ / ^ = 2rE - 2a2t + 0^ . (5.9) 

J^o (/X2 - ai)Y/(/i2 - m5)(^2 - m-e) 

This angle representation is similar to the one found in the case of the defocusing NLS equation. 
(See Alber and Marsden(1994).) 
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Notice that the 0j's are essentiahy the sum of two Kaup type sohtons as t — > ±00. The 
integrals in system (^.^),( |5^ may be evaluated to obtain the following nonlinear algebraic system 
of equations for the 

(■S1625 + •S1526)(S2625 + ^2526) = ^l('Sl625 — Sl526)(s2625 — ^2526) ) (5.10) 

(Sl615 + S1516)(S2615 + S2516) = ^2(si615 " Si5i6)(s2615 " SlSie) , (^-H) 

where Siju = ^/ - mj){ak - mi) and 

Ai = exp[- V(a2 - m5)(a2 - me){2x - 2ait + 9^)] , (5.12) 



A2 = exp[- v^(ai - m5)(ai - mfi){2x - 2a2t + (9^)] . (5.13) 
5.4 Phase Shift Formulas 

When two solitons interact, they normally re-emerge with their initial profile and velocity. However, 
they have shifted ahead or behind where they would have been had there been no interaction at 
all. The amount a soliton shifts is called its phase shift, and the integral equations ( |5.8| ),( |5^ ) can 
be used to compute it. Assume that oi > 02 and define 



M{x) = J^^^. (5.14) 



The phase shifts for the systems (|5.8|) - (|5.9D are calculated by using the following procedure. First, 
consider the reference frame where 9i is a constant, that is 2x — 2ait + 01 = ai for some constant 
ai. Then observe that 62 can be written only as a function of t and ai. Namely, 62 = ai + 2(ai — 
02)^ + ^2 ~ ^1 • Notice that the integrals on the left hand side of (|5.8| ) are exactly the expression for 
single Kaup type solitons. Integrals of this form may be integrated as 

- M(ai) 



log 



,0 2(//2 - ai)^/ (/i2 - rn5){fi2 - me) 2y^{ai - m5)(ai - m^) 



ip + M{ai) 



(5.15) 
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where for clarity we have defined ijP' = M(/i2)^- Then we see that in this frame when t ^ co the 
right side of ( |5.9D goes to infinity so that the left hand side must also grow without bound. On 
evaluation of the integrals we see that this can only happen when — > —M{ai). Similarly when 
t — > — oo, '0 must approach M(ai). Substituting these values for ^2 in (|5.9|) gives asymptotics for 
01 as t ^ 00 and as t ^ —00, respectively. The phase shift for one of the solitons will be the 
difference between the behavior of 61 at minus infinity and its behavior at plus infinity in the frame 
given by 9i = ai. In this way we obtain the shift between phases before and after interaction for 
61 to be 

M{a2) + M{ai) 



log 



M{a2)-M{ai) 



A/(a2 - m^){a2 - ttiq) 
Similarly the phase shift for fi2 along 02 = 02 is computed to be 

M{ai) + M{a2) 



(5.16) 



A2 = — log 

Y^(ai - m5)(ai - me) 



(5.17) 



M{ai) - M{a2) 

These formulas were initially obtained by Matveev and Yavor(1979) by studying asymptotics of 
singular ^-functions. 

5.5 Kink-Anitkink Interaction Solutions 

The interaction of kink and antikink solutions will now be considered. The two-kink solutions are 
constructed by taking the limit of the spectral parameters so that mi m2 ai, m/^ a2, 

and 777-5 ~^ "^6 ~^ 0-3 ■ We examine both the case when the /ij's are both initially on the positive 
branch of the Riemann surface and the case when only one of the //j's is initially on the positive 
branch and the other is on the negative branch. In contrast to the KdV equation, this difference 
in the initial conditions produces qualitatively different solutions. This difference arises because 
the KdV and cKdV equations contain at least one hyperelliptic branch point, while all of the fixed 
points are double points for the two-kink solutions. 
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5.6 Initial values of yUi and fj,2 are chosen on the positive branches of the Rie- 
mann surface. 

In this situation, the angle variables are computed to be 

01 = + 













ilJ'i 


- 02) (mi 


- 03) 














- 02)(A*2 


- 03) 


rtJ-i 




dfii 








- ai)(//i 


- 03) 






dfi2 








- ai){ji2 


- 03) 



+ / . — T^. TT = 2^ - + ^? , (5-18) 

02 = + 

+ / 7 r = 2x-2a2t + e^ (5.19) 

Jul {lJ-2 - ai)[ii2 - a2,) 

These integrals are tractable, and the symmetric polynomials 111IJL2 and fii + /X2 which appear in 
the trace formulas, can be calculated explicitly from the resulting system 

(1 - Ai)jiiiJ,2 + {a^Ai - a2)ini + H2) = Aio^ - a| , (5.20) 
(1 - A2)niH2 + (03^2 - ai){iJ,i + H2) = A2a\ - a\ , (5.21) 



with 



^ (^1 - Q2)(^2 - Q2) ^2(o2-a3)(x-ait+0?) 
(/^i - a3)(M2 - "3) 



(5.22) 



(/^i-a3)(M2 -03) 

These calculations are carried out for the unperturbed Boussinesq equations {K\ = 0) in the next 
section and soliton fusion and fission are discussed. 

5.7 Soliton Fusion and Fission. 

Given this particular deformation of the Riemann surface, solitons can undergo fission or fusion. 
This interesting phenomenon occurs when two separate solitons enter an interaction but only one 
single soliton emerges from this interaction. Since all of the equations under study are invariant 
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under the space-time inversion (x —x,t —t), the reverse process of soHton fission may also 
occur where a single soliton breaks into two distinct solitons at some critical time. Observe that 
soliton fission can be interpreted as an infinite phase shift. This is because the solitons change their 
speed since when they fuse together. Therefore they will be infinitely far from where they would 
have been had there been no interaction. In this sense formulas (Fl6|)-( Fl7|) are still correct since 
in the limit = tuq = 03, M{x) = 1 and so the formulas become singular. 
For this discussion we assume 

that7 = -2^2n+2 

rrii = although the general case is similar. 
Hence the system ( p.20D -( p^ ) can be solved for ^1 + /X2 and is found to be 



TT A^ , \ .01/1+02/2 + 03/3 

= -4(/ii + ^2) = -4 f , f , f > (5-24) 

Jl + /2 + /3 



where 



/i = exp(2aix + 2af t + 9i) , (5.25) 
/2 = exp(2a2X + lajt + 62) , (5.26) 
h = exp(2a3X + 2ait + ^3) , (5-27) 



and 



e\ = ai0? + log[(a2-a3)(^;-ai)(M^-ai)] , (5.28) 

02 = a20? + log[(a3-ai)(^;-a2)(/uO-a2)] , (5.29) 

03 = a30^ + log[(ai-a2)(^?-a3)(M2-a3)] • (5.30) 

The expression ( |5.24| ) coincides with a solution to Burger's equation which describes a confluence 
of shocks (see Whitham(1974)) and has the same form as a solution obtained by using the Hirota 
method. This expression is now analyzed to see why it represents soliton fusion. Suppose oi < 
< a2 < 03 and consider t = 0. We claim that at this instant ( 5.24| ) is a two-tiered kink (sum of 



two kinks), see Figure B^. To see why this is the case, consider 2x < (6*2 — 6*1)7(01 — 02). For 
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these values of x, fi is the largest of the terms /i, /2, and f^. In fact as x — > — oo, the terms 
/2)/3 are negligible in comparison to /i. So for these values of x, /xi + /i2 is nearly constant and is 
approximately oi. Similarly for (6*2 — Oi)/(ai — 02) < 2x < {63 — 62)/ia2 — 0.3), f2 is the dominant 
term and //i + /i2 ~ 02- For the remaining values of x, /s dominates so ^1+^2 ~ 03- 

We now apply this analysis for arbitrary t. In this manner we get that U is essentially constant 



in three regions, Di, D2 and 1)3, of the {x,t) plane as seen in Figure B^. When t is sufficiently 



small, these regions are bounded by the points where the functions /i = /2 and where /2 = /a- At 

^* ^ (Q2 - Q3)(^~2 - Oi) - (qi - «2)(<9""3 - O2) ^^-^ 
2(ai - a2)(a2 - a3)(ai - as) 

the functions /i = /2 = /s for some x*. From this instant on, /2 is never the largest term and so 
for t > t*, the plane is divided into two regions bounded by the points where /i = /s as seen in 



Figure B.7. For more details on this analysis see Whitham's(1974) chapter on Burger's equation. 
As explained above, the regions Di, which contain all the information regarding fission of the 
two-kink solution (times, speeds etc.), are derived from the parameterization of the lines fi = /2, 
/2 = /si and /i = /a- These lines are computed explicitly giving that fi = f2 along the line where 
2x = 2a3t + (6*2 — 6*1)7(01 — 02), /2 = /s along the line where 2x = 2ait + (^3 — 62)/ {0-2 — 0-3) and 
fi = fs along the lines where 2x = 2a2t + (^3 — ^i)/(ai — 03). From this it is possible to predict 
the times at which two solitons experience fission or fusion based solely on the initial phase and 
the values of the three spectrum points Cj. The speeds of the solitons are given by the slopes of 
the lines, and in this way a complete explanation of fission is given. Notice that the arguments are 
quite general and that this approach can be applied directly to the entire class of A'^-component 
systems. 
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5.8 Initial values of /ii and /i2 are chosen on different branches of the Riemann 
surface. 

This case gives a similar problem of inversion as the fusion case, namely 



+ 



Ml 
M? 

2x-2ait + 9^, (5.32) 

I I I rt fin II I I rt f I 'i I 

Ml 



+ 



Ml 







- a2)(/^i 


- as) 




dm 


- 03) 


(/^l 


- ai){fii 
dji2 


- as) 




- ai){n2 


- as) 



- n 7 ^rr T = 2x- 2a2t + . (5.33) 

J^4 (At2 - ai)(/i2 - as) 

These integrals may be evaluated in terms of the log function and give rise to the following algebraic 
equations for the ^j's 

(/ii - ai)(/i2 - as) (^? - ai)(^^ - as) 



(^1 - as)(^2 - ai) (^? - a3)(/i^ - ai) 

(/ii - a2)(^2 - as) ^ (/ig - 02) (/^2 - as) 
(^1 - as)(/i2 - a2) (/i? - a3)(/i^ - 02) 



exp [(ai - as) (2x - 2a2t + 9^)] , (5.34) 
exp[(a2 - as)(2x - 2ait + 6*?)] . (5.35) 



This is a system from which the /ij's can be found explicitly. The solutions are plotted in Figure 



and the corresponding U and W are graphed in Figure |B.10| . By initially choosing different branches 
of this particular Riemann surface we see that the solitons 'change form', i.e. a kink changes 
to an antikink and vice versa. An analytic explanation is given in the next section. This was 
first observed in Alonso and Rues(1992) who noticed this phenomenon asymptotically. Using the 
algebraic geometric construction the finite-time interactions can be analyzed as well. 

5.9 Change of Form of Kinks to Antikinks. 

The solutions obtained from ( |5.34| ) can be put in the following form 

-ai5i - a2ff2 - as5s „„x 

III = ■ ■ , (5.36) 

91+92+ 93 

ai/ll + 02/12 + as/13 rr- or7\ 

= hi+h2 + h, ' 
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where 



91 

92 

93 
hi 
h2 
h3 



exp 
exp 
exp 
exp 
exp 
exp 



-2aix — 2a2a3t — 0162 + log 



-2a2X — laia^t — 02^? + log 



(as - a2)(/ii - ai) 



- "1 



(ai - a3)(/ii - 02) 



/^2 - «2 



-2a3X - 2aia2t - 03(0? + 6^) + log 



(02 - Qi)(Ati - as) 
^0 - aa 



2aix + 20203* + ai^2 + log 
2a2X + 20103* + 026*° + log 



(03 - a2)(/ii - oi) 



A^2 -oi 



(01 - a3)(;U? - 02) 



/x^ -02 



203X + 2oi02t + 03(0? + 9^) + log 



(02 - ai)(/i? - 03) 



f^2 - «3 



(5.38) 
(5.39) 
(5.40) 
(5.41) 
(5.42) 
(5.43) 



Notice the similarity between the form of the solutions ( ^.36 ), ( 5. 37] ) and the solution ( ^.24 ) for U in 
the previous case. Since they have an identical form as that in ( p.24| ), all the analysis from the last 
section applies and we conclude that both fii and /i2 are kinks which experience fission or fusion 
respectively. We find that as t — > —00, fii is an antikink and fi2 decomposes into two kinks. This 
implies that —jU = ^i + ^2 consists of two kinks and one antikink. As t ^ 00, /^i fissions into 
two antikinks and the 2 kinks comprising ^2 fuse into one kink. Therefore /ii + /U2 is the sum of 
one kink and two antikinks. This explains the transformation of kinks to antikinks and vice versa. 



See Figure B.ll to see how these ^ variables combine to form U . Of course the same analysis can 
be performed as above to see this analytically. This is the first time two-soliton solutions of this 
equation have been derived in this simple form. 



6 The SHG Equations 

Solutions of the cKdV may also be transformed into solutions of the SHG equation if n is chosen 
to be —1. This is shown explicitly in Appendix The results from the previous sections can be 
viewed in the context of the SHG equations. First, if k = —1, the initial conditions must be chosen 
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differently or tlie ;Uj's will not be real valued. That is, the cases k = 1 and k = —1 are dual to each 
other in the sense that in the process of the construction of the Riemann surface the cuts in the 
complex plane for the k = 1 case correspond to where real valued solutions lie in the k = —1 case 
and vice versa. This is because the kappa appears under the square root in ( |3.7D so that real values 
of fj, for kappa = 1 correspond exactly to imaginary values of when n = —1 and vice versa. 

Another important difference is that no kink solutions exist for the SHG equations. When 
K, = —1 there is no way to deform the Riemann surface so that all points coalesce piecewise as 
required for kink solitons. Therefore there are always at least two hyperelliptic points left. This 
means that SHG solutions do not have a possibility of either fusing or fissioning and no change of 
form can occur. 

The phase shift formulas for this system are very similar to those in the cKdV hierarchy with 
K = 1, namely 



Ai = ^ log 

y/-{a2 - m^){a2 - me) 

A2 = — ^ log 

V-(ai - m5)(ai - me) 



M(a2) + M(ai) 



M(a2) - M(ai) 
M{ai) + M{a2) 



(6.1) 
(6.2) 



M(ai) - M(a2) 

This is the first time these formulas have been derived. It is remarkable that the formulas are so 
similar to those of the coupled KdV equation since they are not in the same hierarchy of equations - 
the two equations are derived from two different potentials. The phase shift formulas show another 
strength of the algebraic geometric method since the details of the inverse scattering transform 
have not been completed at this time, see Khusnutdinova(1998). u and v are plotted in Figure 
KT% . These funct ions can be transformed into (?i,(?2 of the SHG by transformations in Appendix 



0- 



7 Modified Coupled Dym Equations 
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7.1 Genus-1 solutions for the cDym System 
7.1.1 Periodic solutions. 

A periodic solution of the cDym equation is described by the following differential equation 



n-=i(M-"^i) 



dX , (7.1) 



for particular choice of m^'s. Here X = 2x — 2ct + 9q and integration is carried out on the Riemann 

surface W"^ = — This can be reduced to a standard form by introducing a new variable Y by 
E 



After integration (|7.1|) becomes 



= Y . (7.3) 



Notice that this holomorphic differential is defined on a genus two Riemann surface. To invert this 
integral one has first to consider the following problem of inversion 

Jf^l V/^l<^4(/"l) Jfi^ V/^2C4(/i2) 

■J 1^1 VAtiC4(/"i) J^i°2 V/^2C4(/X2) 
where Ci{^) = nf=i(/^ ~ nij). One needs to rearrange integrals in such a way that to obtain Xi on 
the right hand side of the first integral equation. This yields exact formulas for /xi and ^2 in terms 
of Riemann ^-functions. By fixing ^2 = fn^ and writing fj, = fj,i we resolve the initial problem of 
inversion. (For details see Alber and Fedorov(1999).) 
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7.1.2 Kink solutions. 

Now consider a kink limit by setting mi, m2 ai and m^, — > 02 such that 02 > ai > 0. Integral 
(|7.3|) becomes 

7 ^ = • (^-6) 

(/X - ai)(/x - aa) 

Observe that this is the same problem of inversion as in the cKdV case except that we have Y 

instead of X on the right hand side. After integrating we obtain 

ai(/Uo - 02) - a2(/Uo - ai)exp((ai - a2)y) 

= —f ^ 7 ^ 77 WPT- ■ (7-7) 

(/"o - ^2) - (/"o - oi) exp((ai - a2)Y) 

This gives as a function of Y and as in the cKdV case this is a kink. X is defined in terms of 

Y by ( |7.2| ) once we know /u(y) from ( |7.7| ). The integration in ( [7. 2D may be carried out explicitly. 

Note that dX /dY > and so X(y) is always increasing. By the definition of X it is also clear 

that the range of X(Y) is all real numbers. Therefore the Inverse Function Theorem implies that 

an inverse function Y = Y{X) exists for all values of X and is monotonically decreasing. Therefore 

the graph of ij,{X) = ij,(Y{X)) will have a similar appearance to that of niY), that is, it is also a 

kink. 

After combining numerics for Y{X) with the expression for /u, a description for the kink of the 



cDym system is obtained, see Figure B.13. 



7.1.3 Cusp solution 

For this solution, the limit m\,m2 — > ai and 1713,1714 — > 02 such that 02 < < ai is analyzed. The 
analysis is the same as in the kink case except that now dX /dY changes sign exactly once when 
the branch point Y* is crossed where ^7(Y*) = 0. In this case, Y[X) has two branches begining at 
the hyperelliptic point X* where Y{X*) = Y* . Therefore has two branches and reaches a 



cusp at the point X* . See Figure B.14. 
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7.1.4 Peakon solution 

If the Camassa-Holm shallow water equation is an indication (1994), a peakon may develop in the 
limit mi,m2 ai and 7713,1114 — > 02 with 02 = < ai. The analysis is similar to that of the kink 
case, except now the range of X{Y) is bounded above by some number X*. This means that the 
inverse function Y(X) is defined only for those X < X* . But, it can be defined symmetrically, as 
if the integration were carried out on the negative branch of the square root, and this gives rise to 
a peakon solution. The difference between this and a cusp solution is that in the cusp, 

is infinite at the branch point, while in the Peakon case f^{Y) > for all Y and so at the branch 
point the derivative is finite. 

7.2 Genus 2 solutions 

The algebraic geometric procedure outlined thus far in this paper can be used for other equations 
as well, even when other methods may fail. (For details see Alber and Fedorov (1999).) Using our 
experience with the cKdV system, the case when there are three double points will be considered. 
For the positive branch of the Riemann Surface W"^ = C{E)/E, the problem of inversion can be 
written as 

^^'^^^ +- ,r^^, ^ = dt, (7.9) 



2(mi - ai)(/xi - a2)(Ati - 03)-^^ 2(/i2 - ai){n2 - a2)()U2 - a3)-\//^2 



H 7 TT TT r = dx . (7.10) 

2(mi - ai)(/xi - a2)(Ati - aaj^/JIi 2{fi2 - ai)(|U2 - a2)(/X2 - aajyTx^ 

This inversion problem is very similar to that in the cKdV case, with the exception of the poles 
present in the left hand side of the equation. This case is complicated by the fact that only 
differentials of the third kind appear having simple poles at 01,02,03. Following the procedure 

30 



outlined in Alber and Fedorov(1999), a third variable y is introduced such that 



+ 



dfj.2 



2(/ii - ai)(/xi - a2)(/ii - as)^//!! 2{fX2 - ai)(M2 - a2)(/U2 - az)y/¥2 
Also the normalized differentials of the third kind are introduced 



dy 



(7.11) 



i = 1,2,3 



where = A/a^. Next, consider three points Zi given by 



3 rP- 



i = 1,2,3 . 



The Zi are dependent on x,t,y by 



(7.12) 



(7.13) 



zi = 2^fa{ {x - (a-i + as)* + 0203?/] , 
zi = 2^/a^[x - (ai + 03)* + oiaay] , 
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2V«3 [a:^ - (ai + 02)* + 0102?/] • 



(7.14) 
(7.15) 
(7.16) 



Integrating ( [7.13 ) and putting ^ = ^fji gives 



(6 - «»)(6 - «i)(g3 - 

(6 + ai)(6 + ai)(6 + ai) 



i = 1,2,3 . 



This gives the following system for the symmetric polynomials in ^. 



'^1-e^i -ai(l + e^i) af(l-e^i)^ ^ ^ 



l-e^2 _Q,2(l + e^2) Q,2(i_e^2) 
l-e^3 _Q,3(i + e23) 0,2(1 _e^3) y 



Qf(l + e^i) 
ai(l + e^2) 
y Qi(l + e^3) y 



(7.17) 



(7.18) 



yi-e-^ -asU + e'-'j agU - e-j y \^ 6+6 + 6 / 
/,From this, the expression /Ui + ^2 + /^3 = + ■^i + '?3 be found as a function of zi, 2:2, -^3- 
The determinant of the matrix in the system must be zero to obtain nonzero solutions. This added 
equation gives /i3 in terms of ii\ and /i2 • Then we can find /xi + /i2 as a function of z\,Z2- Then 
this solution must be connected to the {x,i) variables by using ( [7.14 ). 
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Numerics are provided in Figures B.15 and B.16. As expected, the phenomena of change of form 
and fission occurs. A more detailed analysis of fission/fusion for the cDym case will be described 
in a forth coming paper. 
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A Transformations to Related Equations 

In what follows we describe exact connections between solutions of the cKdV and the Kaup equa- 
tions, Boussinesq systems and SHG system. 

A.l Generalized Coupled KdV System 

For this system, assume that k = 1 and r = 0. Observe that if we take n = 1, B takes the form of 
B = E -\-bi. Furthermore, as in ( |B.1[) we see that hi = —u/2 + Ki. Collecting coefficients of order 
1 and respectively in ( |2.9D and substituting the value of 6i shows that u, v satisfy the modified 
cKdV system 



Ut = V - 



\uu' + Klu' 



(A.l) 



1 /// / 

Vt = -gu -vu - 



^uv' + Kiv' . 



(A.2) 



Ki is a small parameter. When Ki = this is the cKdV. 
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A. 2 Generalized Boussinesq System 

To connect u and v to the classical Boussinesq system the following change of variables must be 
made 

u{x,t) = -\U{X,T) , (A.3) 
v{x,t) = ±U{X,Tf-lW{X,T), (A.4) 

where X = x and T = —t/2. Plugging in u and v from ( [A.3| )-( |A^ ) into ( A.lD -( [A^ ) and defining 
7 = —2Ki we see that U and W satisfy 

Ut + Wx + UUx = lUx , (A.5) 
Wt + Uxxx + {WU)x = iWx ■ (A.6) 

Again 7 is a small parameter. When 7 = 0, we have exactly the Boussinesq System. Otherwise 
the system gives rise to a generalization of the Boussinesq system. Such a transformation can be 
found in Sattinger(1995) for example. 

A.3 Generalized Kaup Equations 

To connect U{X,T) and W{X,T) from above to Tr{x,t) and (j){x,t) from the Kaup equations 
(|1.7|) , (|1.8| ) the following change of variables was suggested by Kaup(1972): 

U{X,T) = ^Mx,t) , (A.7) 
W{X,T) = /3-2[l-e7r(x,t)] , (A.8) 

where X = x, T = f3t, f5 = SVS/y/l — 3a, and a = Pj. Then substituting this into ( A.5D -(| 



gives 



VTt = 4>xx + - 3a)d'^(t)xxxx - e{4>x'^)x + OiTTx , (A.9) 

vr = + - a(t)x , (A.IO) 
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which is a perturbed Kaup equations and reduces identicahy to it when a = 0. Summarizing we see 
that every solution of the coupled KdV system yields, using the transformations stated, a solution 
of the Boussinesq and Kaup systems. 

A.4 The SHG System 

For the SHG system we choose k = — 1 in the potential (|2.6D , and define w,^ by the equations 
u = w, and v = i'x/2 + Then the generating equation ( |3.1[ ) becomes 



-2BxE'^ + Bw^E + iBv^x + iBvu^ . (A.ll) 
Now we choose B = b/E, that is we choose m = 1, n = —1, b = 6_2- Then ( [A.ll ) becomes 



Ewt + ^Uxt + ^i^i^t = -^bxxx + 2bxW + ^bxh'x + ^bxV^ 

-2bxE + bwx + i^bvxx + ^bvvx . (A.12) 

Collecting coefficients of E gives 

6=-|, (A.13) 

where we define rj hy rjx = wt- Then collecting coefficients of orders zero and one respectively and 
substituting in the value for b gives 

^i^xt + \vvt = -VxW - \r]Wx (A.14) 

= \'nxxx - ^Vxi^x - Ivxi^^ - jTjVxx - Iw^x ■ (A.15) 

( A.15| ) may be integrated to get 

i'nx)x = iw)'^ + Wx + j Vxirjv - r]x)dx . (A.16) 
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Notice that this equation is satisfied when rjx = r]v. Next we define the new function s by the 
relation vt = — r]w. Then plugging this into ( A.14| ) gives 



rjxW — rjWx H Tjvw = —2rjxW — rjWx ■ 



(A.17) 



Substituting r]x = rjv and canceling like terms yields Sx = sv. These two equations, along with the 
two we defined will determine the SHG system. Summarizing, we have 



r]x = r]u = wt 



rjw 



(A.18) 
(A.19) 
(A.20) 



Next define Q = s ^ and (f)t = r]. This and rjx = wt implies that w = (px- Therefore ( |A.18 ) becomes 



and 



•tlx 



Qx4>t + 



Xt ) 



Sx , Vx 
2^+— ' 

-ivx-v^) = 

s 



(A.21) 



(In Q)xt - (t)x(l>t = - (In s)xt - Tiw , 



Sx\ 

— — nw 
s J t 



-Vt - r]w , 

-{s~^ — rjw) — rjw 



-Q 



(A.22) 



Finally, substitute qi = (v^/2) exp[i((/)/2)] and the real and imaginary parts of the following 
equation correspond to ( |A.21 ) and ( A.22| ), so that 



qixtql - qltQix = -'^{qiQi 



*\2 



(A.23) 
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If we define q2 = —[{lnQ)x + i(t)x] exp(z(/))/4. Then qi and q2 satisfy 

qix = -2q2ql , (A.24) 
q2t = qi , (A.25) 

which is exactly the SHG equation, and qi and q2 are obtained from the fi variables through the 
relation 

cbx = u, (A.26) 

The above transformations were inspired by Khusnutdinova(1998). Summarizing the above 
gives that solutions of the SHG system may be expressed as follows 

qi = - exp ( - I ut ( I ut dx] dx] exp ( - f u dx] , (A. 28) 



^2 = 7 ( -j—^ log ( / utdx ] — iu] exp [ i I u dx ] . (A. 29) 



4 ydx"^ 

Notice that the dependence on v is implicit by the fact that ly = ut [J ut dx) ^ . 

B Trace Formulas 
B.l cKdV System 

In this section, the connection between u, v, and fi for cKdV is derived. This connection is called 



the trace formula for this system. Collect coefficients of in (3T) to see that 2«;6q = from 



which we assume that bo = 1. Now gather coefficients of order n — r + 1 to arrive at 



u 

fei = -^ + i^i, (B.l) 
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for some constant of integration Ki. Next collect coefficients of order 2n — 2r + 1 in (|3.2D to see, 
along with ( |B.1| ), that 



^ 2n+2 

Ki = --^mi. (B.2) 



i=l 

^) then yields the trace formula for u, 

n 2n+2 



^ rui . (B.3) 



u 

i=l i=l 

Now we just need the trace formula for v. For this we collect coefficients of order n — r in ( |3.lD 
to get 

= 2k6^ + 2h[u + hiu' + v' , (B.4) 
which gives upon substitution of hi from ( [B.lD that 

or simply 

= -2k62 + — -K1U + K2, (B.6) 
where -fC2 is some constant and Ki is the same as above. From the definition of B{E) we see that 

l<i<j<n 

All that remains is to find K2. To derive this we collect coefficients of order 2n — 2r in (p.2|). To 
simplify calculations let c„ be the (2n — 2r)*'^ coefficient of the polynomial C{E). Then we see that 

2K{2h2 + hi) + 2u{2hi) + 2v = Cn , (B.8) 
or after substitution of the value of hi that 

4k62 + 2uKi + 2^/^:? + 2v = Cn . (B.9) 

2k 
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Now solving for K2 in ( [B.6D and substituting in v from 

"iv? 

K2 = v + 2Kb2- — + Kiu, (B.IO) 

I - 2Kh2 + ^-uKi- KKf ] + 2Kb2 -^+Kiu, (B.ll) 

= j-'^Kf. (B.12) 

Now from the definition of C{E) we see that 

Cn = 2k mirrij . (B.13) 

l<i<j<2n+2 

B.2 cDym System 

In an analogous manner, we derive the trace formulas for the cDym system. The derivation of the 
trace formula for u is identical to the cKdV case so that 



(B.14) 



n 2n+2 

n = 2' 

i=l i=l 

Only the trace formula for v is different. For this we collect coefficients of order n — r in ( |3.3| ) 
to get 

Ull 

= — ^- + 2h'2 + 2b\u + biu' + v\ (B.15) 
2 

which, upon substitution of 61 from ( |B.1| ) and integrating, gives 

u" 3n^ 

v = -—-2h2 + — -Kiu + K2, (B.16) 

where K2 is some constant and Ki is the same as in cKdV. From the definition of B[E) we see 
that 

h2= ^ Aij/Uj • (B.17) 

l<i<j<n 
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All that remains is to find i^2- To derive this we collect coefficients of order 2n — 2r — 1 in ( p.4[ ). 
To simplify calculations let c„_i be the (2n — 2r — 1)*** coefficient of the polynomial C{E). Then 
we see that 



-b'l + 2(262 + bl) + 2n(26i) + 2v = c^-i , 



(B.18) 



or after substitution of the value of bi that 



u" 2,v? 

y + 462 - — + 2uKi + 2Kf + 2V = Cn-l 



(B.19) 



Now solving for K2 in ( [B.16| ) and substituting in v from ( B.19 ) gives 



3u^ u" 
K2 = v + 2b2 - — + Kiu + — , 

= (|-262 + ^-ni^i-i^?-^) + 262 



iu' u" 
— + K,u+-, 



Kl. 



(B.20) 
(B.21) 
(B.22) 



^Prom the definition of C{E) we see that 



Cn-l = 2 ^ 1 
l<i<3<2n+2 



(B.23) 
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Figures 




Figure B.l: Periodic /i solving ( |4.7| ) at t = in the genus-1 case. This corresponds to distinct roots 
rrii in (p.5|). This plot uses mi = 0, m2 = 3, = 2, = 1. From this function fi, solutions u,v 
of the cKdV equation may be obtained by ( 3.12| ),( prT3 ). Then U, W from the Boussinesq equation 



may be found using (A. 2) and 



44 



40 




35 






30 








f 


M 1/1/ 


25 






20 






15 






10 






5 
























5 ' ' 1 ' ' ' 1 ' ' 1 1 

10 8642 02468 10 



Figure B.2: Genus-one solutions U and W to the Boussinesq equation obtained from ( 4.15 ) using 

a = —2, TTT-s = 1, 7714 = 3. 




Figure B.3: U and W shown here are genus-one solutions to the Boussinesq equation for two double 
roots, t = 0, ^0 = 0, c = —4. This graph highlights the simple relationship between U and W 
which will be important in understanding two soliton interaction. 
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Figure B.4: Genus-two solutions U and W to the Boussinesq equation for one double root derived 



from (|5.5D , (|5.6|) and the trace formulas. This corresponds to a Kaup type soliton on a quasi-periodic 
background. Here t = 0, a = — 5,7713=0,7714 = 1,7775 = 3,7776 = 3.5, ^i(O) = —4.99, /i2(0) = 2. 




Figure B.5: These graphs show fii and fi2, genus-two solutions with two double roots, for al = —5, 
a2 = —14,7775 = 5,7776 = 7.5,t = —.07. This illustrates two Kaup type solitons. 
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Figure B.6: This figure shows genus-two solutions U and W to the Boussinesq equation derived 
from Q, Q and (jA^, (Oh. Here al = -5, a2 = -U,m5 = 5,m6 = 7.5,t = -.07. . . .07. 




Figure B.7: Areas Di, D2, where away from the boundaries, U from the Boussinesq equation 
has a constant height of 01,02,03 respectively. Thus for a time t < t* , U is the sum of two kinks 
while for t > t*, it is only one kink. This illustrates soliton fusion. 
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Figure B.8: Genus-two solutions U and W to the Boussinesq equation with two double roots. Here 
al = —2, a2 = 2,03 = 5,t = — A. . . .4. Both fi variables are initially chosen on the positive branch 
of the Riemann surface. Observe that each variable U and W undergoes soliton fission. 




Figure B.9: Genus-two /ii and fi2 from ( 5.36 ) and ( ^.37 ) illustrating when fii and fi2 are initially 
on different branches of the Riemann surface. The parameters chosen for these graphs are al = —2, 
a2 = 0,a3 = 2,t = —2 ... 2. Notice how similar these graphs are to the fission/fusion of U in figure 
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Figure B.IO: Genus-two U and W from the Boussinesq equation derived from ( 5.36D , (|5.37|) and 
([All) , ([A^ ). Here al = -2, a2 = 2,a3 = 5,t = -A... .4. Here the variables from (|5.36D and 



( 5.37| ) are chosen to be on the opposite branches of the Riemann surface. Notice the change of form 
in U as kinks become antikinks and vice versa. 




Figure B.ll: A contour plot of a genus- two solution U of the Boussinesq equation for the parameters 
al = —1, a2 = l,a3 = 2,t = —3. . .3. This shows how the /j, variables combine to form different 
parts of the same kink. 
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Figure B.12: A plot of u and v, genus- two solutions derived from the variables in (|5.8|) and (^.91) , 
except that k = —1, using the trace formulas. These can be transformed into the SHG equation 
using transformations in appendix A. Here (k = —1) for al = —1, a2 = l,m5 = 2,mQ = —2,t = 0, 
and the variables are initially chosen on opposite branches of the Riemann surface. 



Figure B.13: A genus-1 kink solution, u{x,0), for the coupled Dym equation. The parameters used 
were two double roots al = 1, a2 = 2, and t = 
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Figure B.14: A genus-1 cusp solution for the cDym equation, u for two double roots al = — 1, 
a2 = 1, and t = 0. 
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Figure B.15: Genus-two solution u for the cDym equation with three double roots al = .1, a2 = 
.5,03 = 1,* = — 12...12. Here the /x variables from which u was derived were initially chosen to lie 
on the same branch of the Riemann surface. In this situation soliton fission occurs as in the ckdv 
case. 
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Figure B.16: Change of form for the Coupled Dym equation is seen for the genus-two solution u 
with three double roots al = 1, a2 = 6,03 = 10,i = —.3 ... .3. The variables from which u was 
derived were initially chosen on opposite branches of the Riemann surface. 
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